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Shape Parametrisation 
This invention relates to a method of and apparatus for shape parametrisation of digital 

data. 

US patent specification No 3069654 describes a method of finding the parameters of 
straight lines in an image. Each point in the image space is mapped into a gradient, intercept 
5 (m,e) parametric transform space to produce lines representing all possible gradients and 
intercepts of lines passing through that point. Thus, a point in the image space is 

mapped into a line satisfying the equation c = y, - mx,- in the parametric transform space- A 
maximum, or intersection, in the transform space is detected and determined to represent a 
line in the image space, since an intersection at (m/.c/) of a plurality of lines in the transform 

10 space denotes a corresponding plurality of colinear points in the image space all lying on the 
line having the equation y = m/x + cj. However, when the image space contains a plurality of 
lines, additional spurious maxima or intersections are produced in the transform space and are 
erroneously determined to represent lines in the image space. Furthermore, using the method 
described above, it is not possible to determine whether or not a maximum in the transform 

15 space represents a single line segment or a plurality of disconnected colinear points in the 
image space. A further practical problem with the method described above is that the m and 
c parameters are unbounded and therefore some groups of points in the image space cannot 
be represented in a bounded transform space. 

In view of the last-mentioned disadvantage of the method described above, it is known 
20 to map points in the image space into an angle, radius (0, r) normal parametrisation transform 
space or "sinogram*, rather than an (m^c) transform space. Each point in the image space 
is mapped into a sine curve in the sinogram representing r and 0 for all possible lines in the 
image space passing through the point, where r is the algebraic distance from the origin to 
the line along a normal to the line and 6 is the angle of the normal to the x-axis. Typically, 
25 the limits of B are -tc < 9 < w, in which case for a bounded square image space of size LxL, 
the limits of r are -y/2L <r< y/lL. However, this development of the first-mentioned method 
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still suffers from the problems of spurious maxima and disconnected colinearities. 

The prior art may be classified as providing a method of shape parametrisation of digital 
image data comprising the steps of: 

3 0 transforming the image data into a parametric transform space; and 

extracting shape characterising parameters from the transform space indicative of a 
shape in the image represented by the image data. 

The present invention seeks to overcome the problems of spurious maxima and discon- 
nected colinearities associated with prior art. 

35 The method provided by a first aspect of the present invention is characterised in that 

the extraction step includes the step of detecting at least one particular shape indicative 
distribution of data in the transform space. 

Due to the problems associated with the prior art, it was virtually impossible to provide 
automatic shape parametrisation of the image data even when the data represented very simple 
40 shapes, and prior knowledge of the shape was needed, together with human interface to select 
those of the shape characterising parameters representing a "proper" shape . However, the 
first aspect of present invention provides a more reliable method which can be performed fully 
automatically for more complex shapes than has been possible using the prior art methods. 

Preferred features of the method and other aspects of the invention are set out in the 
45 claims. 

There follows a description by way of example of specific embodiments of the present 
invention, reference being made to the accompanying drawings, in which: 

Figure 1 is a schematic diagram of an apparatus of one embodiment of the invention; 

Figure 2 is a perspective view of an object, the processing of an image of which is 
50 described below; 

Figures 3 to 5 are representations of the image of the object after various processing 
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operations; 

Figures 6 and 7 illustrate the mapping of a single point and three colinear points, 
respectively, from an image space to a sinogram; 

55 Figures 8 and 9 are graphical representations of distributions of curves in the sinogram; 

Figure 10 is a matrix of mask values used in detecting distributions of data in the 
sinogram indicative of a line segment in the image; 

Figures 11 and 12 are digital representations corresponding to Figures 8 and 9, respec- 
tively; 

60 Figures 13 and 14 illustrate the mapping of a circle from an image space to a sinogram; 

Figure 15 is a graphical representation of data intensity across a belt produced in the 
sinogram at the location indicated by the lines XV - XV in Figure 14; 

Figure 16 is a digital representation of the data shown in Figure 15; 

Figures 17 and 19 are mask values for use in detecting data in the sinogram representa- 
6 5 tive of a circle; 

Figures 18 and 20 show the data of Figure 16 after convolution using the masks of 
Figures 17 and 19, respectively; and 

Figures 21 and 22 show the data of Figures 18 and 20, respectively after further pro- 
cessing. 

70 Referring to Figure 1 of the drawings, a camera 10 outputs a digitised video signal 

representing the object illustrated in Figure 2, which is stored as a frame in a frame store 12. 
Typically, the frame size is 256 pixels by 256 pixels, and each pixel has eight bits and so can 
store values in the range 0 to 255. A parallel processor 16, such as a linear array processor 
as described in U.K. patent specification No. 2129545B, then performs an edge detection 

75 operation on the stored frame using a Sobel-type operator in a known manner to produce 
an image as represented in Figure 3 which is stored as a frame in the frame store 12. The 
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image is then subjected to a thresholding operation by the parallel processor 16 to produce a 
binarised image in which each pixel either has a predetermined low value, for example zero, 
or a predetermined high value, for example 255. The binarised image, which is represented 
80 in Figure 4, is stored in the frame store 12. The edges of the image are then thinned, and 
isolated points are removed by the parallel processor 16, and the resulting binarised image, as 
represented in Figure 5, is stored in the frame store 12. 

Once the binarised image has been formed, the edge points are mapped from the image 
space into a sinogram or angle, radius normal parametrisation space (0, r) by a host computer 
8 5 20, and the sinogram is stored as a frame in the frame store 12. Referring to Figure 6, each edge 
point P t at coordinates (x,-, y, ) in the image space is transformed to a sine curve C, representing 
the angles and radii (0, r) of the normals to all possible lines passing through the point (x,, y,) 
in the image space. Thus, the sine curve d satisfies the equation 

r=(x?+j/r)*cos{^tan- 1 (-)}- 

90 By way of example, lines f x , U are shown in the image space of Figure 6 which produce points 
at n) and (02, r 2 ) in the sinogram. 

Figure 7 shows how three points Pi, ft, ft in tbe image space are transformed into three 
sine curves C lf C 3 ,C 3 in the sinogram. Since the three points P u P 2t P^ are colinear, the sine 
curves C^Gj, C 3 intersect at a single point / in the sinogram. The coordinates [6 L ,r L ] of the 
95 intersection / in the sinogram give the angle and length of the normal in the image space 
which define the line L on which the three points ft, ft, ft lie. The line satisfies the equation 
y = cot[—$i)x + r^cosecB. 

All of the points in the image space are transformed into the sinogram in the manner 
described above, and the sinogram is stored digitally as a frame in the frame store 12. In 
100 the case of the three points shown in figure 7, the pixel in the sinogram corresponding to 
the intersection / would have a value of 3. Obviously, in practice, many more points can be 
processed, producing many curves in the sinogram and higher pixel values than occur in the 
simple example described above. 
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In the prior art, mere maxima in the sinogram are detected. However, in accordance 
105 w * t ' 1 tnis embodiment of the invention, more specific distributions of data in the sinogram 
are sought. It has been noted that curves in the sinogram representing a continuous line in 
the image space resemble, at and around the intersection or maximum in the sinogram, a 
butterfly shape with the wings of the butterfly extending in the $ direction, as represented 
graphically in Figure 8 and digitally in Figure 11, whereas a discontinuous line, although 
110 producing a maximum in the sinogram, has a less dense packing of the curves which form 
the wings of the butterfly, as represented graphically in Figure 9 and digitally in Figure 12. 
Thus, by applying an appropriate mask to groups of pixel data in the sinogram, it is possible 
to discriminate between a maximum having a butterfly shaped distribution surrounding it 
and other maxima, and therefore it is possible to detect points which, in the binarised image, 
115 represent a continuous line. 

Figure 10 shows a 3 pixel x 3 pixel mask for detecting butterfly-shaped distributions 
in the sinogram. The sinogram in the frame store 12 is convolved with the mask using the 
parallel processor 16, and the results of the convolution operation are stored as a frame in 
the frame store 12. More specifically, considering the upper three rows of pixels in a frame, 

120 if the pixel values /ru» to /„. 25 5, A.o to I lm2 s& and J 3 .o to /3.25s, the mask is applied to each 
3x3 group of pixels, that is / 0 .,-x, /o.,*, WiJ 'i..-i> A.i, '2.1-1, h.i and where 

1 < t < 254, and the pixel values are multiplied by the corresponding mask values; the products 
are summed. Preferably, the sum is then divided by (1+ the value of the middle pixel of 
the group), and then the result is used as a pixel value J lmi at a corresponding location in 

125 a further frame store 12. Thus, for the mask shown in Figure 10, the pixel value Ju = 
{(/<u x -2) + (J u -i x 1) + (Ji., x 2) + (J,.,- +l x 1) + (J 2J * x -2))/{I u + 1}. The mask is applied 
to 254 3 x 3 groups of pixels in the first three rows of the frame section 18 to produce 254 
results. The operation is then repeated for the second to fourth rows of the frame section, for 
the third to fifth rows, and so on, finishing with the 254th to 256th rows. 

130 Referring to the examples of pixel data shown in Figures 11 and 12, the Figure 11 

group gives a result after convolution of {(1 x 0) + (0 x -2) + (1 x 0) + (7 x 1) + (7 x 2) + (7 x 



1) + (0 x 0) + (0 x -2) + (2 x 0)}/{7 + 1} = 3 whereas the Figure 12 group gives a result of 
{(3xO) + (lx-2) + (3xO) + (4x 1) + (7 x 2) + (3 x 1) + (4 x 0) -f (2 x -2) + (3 x 0)}/{7+1} = 1 

It will be noted that a 3x3 group of pixels each having the same value will produce a 
135 result of zero after convolution. 

After the convolution operation, the host computer 20, selects those of the convolution 
results above a predetermined threshold value, say 2, and representative of a continuous line 
and, from the location of the pixel in the frame section 18, determines the parameters of the 
line represented by that pixel. The convolution result of 3 for the group of pixel data shown 
140 i n Figure 11 produces an indication of a continuous line, as compared with the result of 1 for 
the group of Figure 12, which is not treated as indicating a continuous line, despite the centre 
pixel of that group being a maximum. 

Starting with the detected line having the highest convolution result, the computer 20 
then compares the detected lines with the original image stored in the frame store 12 and 
145 determines the locations of the end points of the detected lines. 

Thus, the r, 0 or m, c and end point parameters of all the continuous lines in the image 
are determined. 

Whilst the above description has been confined to the parametrisation of straight lines, 
it is also possible to determine the parameters of circles, part-circles and other conic sections 
150 in the image. 

Referring to Figure 13, the points on a circle C in the image space are mapped in the 
same way as described above and produce in the sinogram shown in Figure 14 a belt 26 of 
sine curves. For an even distribution of points lying on the circle in the image space, the 
distribution of sine curves across the belt in the sinogram is not even, but rather the curves 
155 have maximum intensities at the edges 28, 30 of the belt and the intensity exhibits an inverse 
square root dependence near the edges. Figure 15 is a plot, by way of example, of the intensity 
/ across the belt 26, and figure 16 is a digital representation of the intensity. 



The parameters of the circle in the image space can be determined from the size and 
phase of the belt in the sinogram. Specifically, the radius p of the circle is equal to one half of 
160 the width of the belt in the r direction and the coordinates (x.„sa.) of the centre of the circle 
are equal to the values of r at the centre of the belt at 6 = 0 and 6 = tt/2, respectively. In 
Figures 13 and 14, the scale of the r axis is one half of that for the x or y axes and so 2p in 
Figure 14 appears to be the same distance as p in Figure 13. 

In the case of an arc of a circle in the image space, the centre (x,., t/o) and radius p of 
165 curvature can be determined from the values of r and 9 at three locations on the edges of the 
belt by forming three simultaneous equations from the following two equations: 

r = xn cos 8 + yo sin 6 + p U) 

r = X( t cos 0 + j/n sin 6 - p ( 2 ) 
where equations 1 and 2 are used for values of r and $ on the upper and lower edges, respectively, 
170 of the belt. 

Two 1 x 3 masks are used to detect the edges of belts, the mask of Figure 17 having 
values (2,-1,-1) being used to detect lower edges, and the mask of Figure 19 having values 
(-1,-1,2) being used to detect upper edges. In the convolution process the upper and lower 
edge masks are traversed in the -r direction of the sinogram along each column of pixels with 

175 a step of one pixel, and at each stage the group of three pixel values in the respective column 
of the sinogram are multiplied by the corresponding mask values and summed and placed in a 
location corresponding to the middle pixel of the group in a further frame. The part-column 
of pixel values shown in Figure 16 after convolution using the lower edge mask are changed 
to the values shown in Figure 18, and after convolution using the upper edge mask take the 

180 values shown in Figure 20. By way of further explanation, referring to Figures 16, 19 and 20, 
the first three pixels' values (0,0,0) on the left in Figure 16 are multiplied by the corresponding 
values (- 1, -1, 2) of the mask and summed to produce a result of (0 x - 1) + (0 x - 1) + (0 x 2) = 0 
which is stored as the value of second pixel, as shown in Figure 20. The process is repeated 
for the second to fourth pixel values of Figure 16 and the result of zero is stored as the value 
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185 of the third pixel, as shown in Figure 20. For the third to fifth pixel values (0,0, 16) in Figure 
16, the result is (0 x -1) + (0 x -l) + (16 x 2) = 32, which is stored as the fourth pixel value as 
shown in Figure 20. The process is repeated all the way down the pixel column, and similar 
processes are carried out by the parallel processor 16 simultaneously for all the pixel columns. 

It will be noted that negative results of the convolution process, such as obtained when 
190 the fifth to seventh pixel values of Figure 16 are operated on by the upper edge mask of Figure 
19, are stored as zero. 

Detection of belts in the sinogram may be carried out by comparing the convolution 
results obtained by using the masks of Figures 17 and 19. Comparing the example pixel values 
in Figures 18 and 20, it can be seen that an upper edge of the belt is indicated by a moderately 
195 high value (18) in the Figure 18 and a higher value (32) at a generally corresponding location 
in Figure 20. Similarly, a lower edge of the belt is indicated by a high value (32) in Figure 18 
and a moderately high value (18) at a generally corresponding location in Figure 20, 

Alternatively, detection can be carried out after further processing of the values in 
Figures 18 and 20 by dividing each value by (1 + the corresponding value in Figure 16) to 
200 produce the values shown in Figures 21 and 22. It will be seen that in Figure 21, the lower 
edge of the belt is indicated by a high value, and the other values are all low, whereas in Figure 
22, the upper edge of the belt is indicated by a high value amongst low values. 

The locations of the pixels indicating the edges of the belt are used by the computer 
20 in determining the sizes and phases of any belts in the sinogram indicative of circles, the 
205 parameters of which can then be determined in the manner described above with reference to 
Figures 13 and 14. 

Ellipses in the image space can be detected in a similar manner to circles. However in 
the case of an ellipse, the belt produced in the sinogram has a varying width. The length and 
width of the ellipse along the major and minor axes are equal to the maximum and minimum 
210 widths, respectively, of the belt in the r direction. The coordinates (sn, of the centre of the 
ellipse are equal to the values of r at the centre of the belt 0 = 0 and 9 = respectively. The 
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tilt of the ellipse, that is the angle of the major axis to the x-axis, is equal to the value of B at 
the maximum width of the belt. 

Other conic sections produce belts or distinctive distributions of curves in the sinogram 
215 which can be detected using an appropriate mask and convolution process. 

Once the parameters of lines and curves in the image space have been determined, they 
can be stored away in a library using much less memory than would be required to store a 
whole frame of image data. Thus the method described above can be used to read, encode 
and store automatically and efficiently 2-D images such as technical or architectural drawings 
220 or electrical circuit diagrams. The method can also be used for example in interferometry to 
detect and provide the parameters of interference fringes. In robotics, the two dimensional 
spatial relationships between lines and curves can be determined by the computer and then 
compared with stored 2-D data. 

In a development of the method described above, two cameras are situated to provide 
225 a pair of stereoscopic images which are both processed as described above in parallel with 
each other. The two sets of image data are then matched by the computer so that a three 
dimensional representation of an object viewed by the cameras can be determined. It is then 
possible to process the 3-D representation and compare it with stored 3-D data to determine 
the identity and location of the object, thus providing automatic robotic vision. 

230 By using a parallel processor as described above, it is possible to process the image data 

one line at a time and thus quickly. If Wgh-speed processing is not required, however, serial 
processing of the image data may be performed. 
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Claims 

1. A method of parametrisation of shapes in images represented by image data, com- 
235 prising the steps of: 

transforming the image data into a parametric transform space; and 

extracting parameters from the transform space indicative of a shape in the image 
represented by the image data; 

characterised in that: 

240 the extraction step includes the step of detecting at least one particular shape indicative 

distribution of data in the transform space. 

2. A method as claimed in claim 1, wherein the transformation step produces bounded 
parameters in the transform space. 

3. A method as claimed in claim 1 or 2, wherein the transform space into which the 
245 image data is transformed is a sinogram. 

4. A method as claimed in any preceding claim, wherein the extraction step includes 
the step of convolving the transform space with at least one mask. 

5. A method as claimed in any preceding claim, wherein the or one of the particular 
distributions of data which is detected is a distribution of data around a point in the transform 

250 space indicative of a line in the image represented by the image data. 

6. A method as claimed in claim 5, when appendant directly or indirectly to claim 3, 
wherein the line indicative distribution of data which is detected is a distribution resembling 
a butterfly shape substantially of the type described in the description. 

7. A method as claimed in claim 6 when appendant directly or indirectly to claim 4, 
255 wherein the mask operates to add to the value of the data at each point in the sinogram at 

least a portion of the values of the data to either side of that point in the angle direction of 
the sinogram and to subtract from that value at least a portion of the values of the data to 



either side of that point in the radius direction of the sinogram. 

8. A method as claimed in claim 7, further comprising the step of dividing the value 
260 0 f tne data at each point after operation of the mask with a value related to the value of the 

data at that point before operation of the mask. 

9. A method as claimed in any of claims 6 to 8, wherein the extracting step includes 
the step of extracting the position in the sinogram of the centre of the butterfly shape. 

10. A method as claimed in claim 9, wherein the extracting step further includes the 
265 step of determining parameters of a line from the extracted position in the sinogram. 

11. A method as claimed in any preceding claim, wherein the or one of the particular 
distributions of data which is detected is a distribution of data extending across the transform 
space indicative of a curve in the image represented by the image data. 

12. A method as claimed in claim 11, when appendant directly or indirectly to claim 3, 
270 wherein the curve indicative distribution of data which is detected is a distribution resembling 

a ridge at an edge of a belt. 

13. A method as claimed in claim 12, wherein the curve indicative distribution of data 
which is detected is a distribution resembling a pair of ridges at opposite edges of a belt. 

14. A method as claimed in claim 13, when appendant directly or indirectly to claim 4, 
275 wherein different such masks are used for detecting the two edges of the belt. 

15. A method as claimed in claim 14, wherein one of the masks operates to subtract 
from the value of the data at each point in the sinogram at least a portion of the value of 
the data to one side of that point in the radius direction of the sinogram and the other mask 
operates to subtract from the value of the data at each point in the sinogram at least a portion 

280 of the value of the data to the other side of that point in the radius direction of the sinogram. 

16. A method as claimed in claim 15, further comprising the step of dividing the value 
of the data at each point in the sinogram after operation of the masks with a value related to 
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the value of the data at that point before operation of the masks. 

17. A method as claimed in claim 15 or 16, and further comprising the step of comparing 
285 w i fc h each other the results obtained by operation of the two masks. 

18. A method as claimed in any of the claims 12 to 17, wherein the extracting step 
includes the step of extracting the positions in the sinogram of at least three points on the 
edge or edges of the belt. 

19. A method as claimed in claim 18, wherein the extracting step further includes the 
290 step of determining parameters of a circle or arc of a circle from the extracted positions in the 

sinogram. 

20. A method as claimed in claim 18 or 19, further comprising the step of determining 
whether the width of the belt in the radius direction of the sinogram is substantially constant. 

21. A method as claimed in any of the claims of 12 to 17, wherein the extracting step 
295 includes the step of extracting the positions in the sinogram of at least five points on the edge 

or edges of the belt. 

22. A method as claimed in claim 21, wherein the extracting step includes the step of 
determining parameters of a conic section or part-conic section from the extracted positions 
in the sinogram. 

300 23. A method as claimed in claim 22, further comprising the step of determining the 

radius parameters of the centres of the belt in the sinogram at angles of nn and (n-f (where 
n is an integer), the maximum and minimum widths of the belt in the radius direction of the 
sinogram, and the phase of variations of the width of the belt. 

24. A method as claimed in any of claims 7 to 10 and any of claims 14 to 17, wherein 
305 butterfly-shaped distributions are not detectable by operation of the masks for detecting edges 
of belts, and the belt-like distributions are not detectable by operation of the mask for detecting 
butterfly-shaped distributions. 
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25. A method as claimed in any preceding claim, and further comprising the step of 
comparing data representing the shape characterised by the extracted parameters with the 

310 image data to determine further parameters of the shape. 

26. A method as claimed in any preceding claim, wherein the method is performed on 
two stereoscopically related sets of image data, and further comprising the step of matching 
the shape indicative parameters extracted for each set of image data. 

27. A method as claimed in any preceding claim, and further comprising the step or steps 
315 of performing the edge detection operation and/or a binarising operation and/or a thinning 

operation on the image data prior to the transformation step. 

28. A method as claimed in any preceding claim, wherein the image data is provided 
by a digital image signal. 

29. A method of shape parametrisation of shapes in images represented by image data, 
320 substantially as described in the description with reference to the drawings. 

30. An apparatus specially adapted to perform the method of any preceding claim. 

31. An apparatus for parametrisation of shapes in images represented by image data, 
comprising: 

means to receive image data; 

325 means to transform the image data to parametric data; 

means to detect at least one particular shape indicative distribution of data in the 
parametric data; and 

means to output parameters related to such a detected distribution. 

32. An apparatus as claimed in claim 31, and including means to store the parametric 
330 data in the form of a sinogram. 

33. An apparatus as claimed in claim 32, wherein the detection means includes means 



defining at least one mask and means to convolve the sinogram with the mask. 

34. An apparatus as claimed in claim 33, wherein the mask is operable to add to the 
value of data at each point in the sinogram at least a portion of the values of the data to either 

335 side of that point in the angle direction of the sinogram and to subtract from that value at 
least a portion of the values of the data to either side of that point in the radius direction of 
the sinogram. 

35. An apparatus as claimed in claim 34, wherein the detection means is operable to 
divide the value of the data at each point after operation of the mask with a value related to 

340 the value of the data at that point before operation of the mask. 

36. Apparatus as claimed in any of claims 33 to 35, wherein one such mask is operable 
to subtract from the value of the data at each point in the sinogram at least a portion of the 
value of the data to one side of that point in the radius direction of the sinogram and another 
such mask is operable to subtract from the value of the data at each point in the sinogram at 

345 least a portion of the value of the data to the other side of that point in the radius direction 
of the sinogram. 

37. An apparatus as claimed in. claim 36, wherein the detecting means is operable to 
divide the value of the data at each point in the sinogram after operation of the masks with 
a value related to the value of the data at that point before operation of the masks. 

350 38. An apparatus as claimed in claim 36 or 37, and wherein the detecting means is 

operable to compare with each other the results obtained by operation of the two masks. 

39. An apparatus as claimed in any of claims 31 to 38, and further comprising means 
to store the image data, the detecting means being operable to compare data representing 
a shape characterised by the parameters related to the detected distribution with the stored 

355 image data to determine further parameters of the shape. 

40. An apparatus as claimed in any of claims 31 to 39, comprising a further such data 
receiving means, transforming means, detection means and output means so that the apparatus 
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can parametrise two stereoscopically related sets of image data, and further comprising means 
for matching the output parameters for each set of image data. 

360 41. An apparatus as claimed in claim 40, further comprising a pair stereoscopically 

arranged cameras for feeding image data to the data receiving means. 

42. An apparatus for shape parametrisation of shapes in images represented by image 
data, substantially as described in the description with reference to the drawings. 
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